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Abstract: We propose a two-dimensional exciton model showing the possibility of two-exciton 
pairing in cuprates, both in singlet and in triplet states. The former is assumed to be responsible 
for the superconductivity and the latter for the unusual normal state properties. The ground 
state is found to be a projection of the Bose condensation of N/2 mobile bosons with a 
Tj projection operator which is a generalization of the Gutzwiller projection operator. The phase 

diagram of the model shows insulator-superconductivity and superconductivity-metal 
transitions with doping. 
PACS No: 74.20.Mn, 74.25.-q 

It has been more than one decade since the discovery of high-Tc cuprate superconductors. 
The core mechanism of Cooper pairing in those materials that is essential to 
superconductivity is still unclear. The coherent length for the cuprate superconductors is 
extremely short (usually about <^~14A), orders smaller than that of conventional 
superconductors (| ~ 10000 A). The Coulomb repulsion in such a short coherent length is 

very strong, so it is unlikely to have an attractive interaction between two holes. However, it 
was realized that an attractive interaction is not the only possible way to the formation of 
Cooper pairs. Anderson [1] first showed that the resonating valence-bond (RVB) states for a 
strongly correlated system are related to an appropriate projection of the BCS function. It is 
believed that the Hubbard model or t-J model should be enough to describe the high-Tc 
superconductivity [2]. The exchange energy J=t 2 /U~0.4ev (t is hoping energy and U is on-site 
Coulomb interaction) required to break a valence-bonded pair, however, leads to a very high 
mean-field temperature T M f~1000K»Tc. It indicates at least that those models may be 
insufficient to describe the low energy physics much below the T M f. Recently, Zhang [3] 
proposed a phenomenological theory about the SO(5) symmetry of superconductivity and 
antiferromagnetism. It is actually a generalization of the SO(4) symmetry for the Hubbard 
model [4]. New operators k and 7t + symmetry operators carrying charge 2 and spin 1 are 
introduced, relating to triplet states by definition. An extended Hubbard model with the 
SO(5) symmetry can be constructed when the complicated infinite long-range interaction 
terms are introduced [5]. It is plausible that the triplet states and the long-range Coulomb 
correlation might be relevant to the superconductivity in cuprates. In this paper we propose 
a two-dimensional exciton model that includes explicitly the Coulomb interaction among 
acceptors and doping holes. The assumption of two-dimensionality bases on the experiments 
that show the incoherent transport of holes in the z-axis direction in the normal state [6] and 
even in the superconducting state [7]. 
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The most important structure element of the high-Tc superconductors (both 214 compounds 
and 123 compounds) is the Cu-O planes. The metal or rare earth ions are sandwiched by 
those planes. Doping generates holes in the CuO planes and acceptors in the inter-spacing of 
two CuO planes. Both quantum chemistry calculations and experiments find that the doping 
holes have 85% O 2p and 15% Cu 3d character. Here we distinguish the doping holes and 
the d-holes 3d 9 for the Cu 2+ ions. The latter exist even without any doping (hereafter holes 
refer to the doping holes). We use a continuous representation of Hamiltonian instead of the 
tight-binding (TB) models in advantage of taking account of large-range Coulomb 
correlation. The Hamiltonian for a system of doping holes may be written as 

where iV is the number of holes and M is the number of acceptors (N < M). The first term 

on the right side is kinetic energy corresponding to the hoping term in the TB models, m* is 
the effective mass of holes (which may depend on momentum as well). The second term is 

the interaction between holes and acceptors V,(r) = -1/ ^Jd 2 + (r-R ; ) 2 ? where d is the 
distance from acceptor to CuO plane, R ; is the position in CuO plane projected from 
acceptor / . The third term is the Coulomb interaction between holes, in which e is an 
effective dielectric constant in CuO plane. The d-holes are not included in the Hamiltonian 
above, but they are indirectly relating to the effective mass and the effective dielectric 
constant. Notice that those holes are confined in the two-dimensional CuO plane. 

It is easy to obtain the solution for a one-hole system of Hamiltonian H(l,l) . It is known that 

in a two-dimensional space there is always a bound state in a potential well no matter how 
shallow it is[8]. Therefore, this hole is in a bound state. A numerical solution is carried out to 
solve the problem. We choose a linear variational approach [9]. The trial wave function may 

K 

be written as 4> = ^c,0, where {(/) , } are a set of basis functions. The energy expectation 

i 

\&H ®dz 

value is E = ^—z — , where H Q is /f (1,1) . Both the eigen wave function and the eigen 

<&*<&dt 

value of energy may be obtained from equation det(H - £S) = 0, where H and S are 
respectively the KxK matrix of the one-hole Hamiltonian and the overlap of basis 
functions. In our calculation, a basis set with three Gaussian functions 
0j =exp(-a r 2 ), z = 1,3 is chosen, and the distance d is set as 4.0 atomic unit (au.). The 
exponents a, may be optimized to give the lowest eigen value E . The results are listed in 
Table 1, in which we include also the solution for a hydrogen atom under the same approach 
for comparison. The coefficient c 3 is two orders smaller than c, and c 2 , indicating a good 
convergence achieved by three Gaussian basis functions. The expectation value of radius 
< r > for the hole is very large due to the wide half-width of potential well. We may call the 
bound state of the hole as a planar exciton. 
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The two-hole system of Hamiltonian H(N = 2, M = 2) may be solved with the help of the 
solution for the one-hole system. The motion of two holes is correlated through the long- 
range Coulomb interaction l/er 12 . This system has two spin states: singlet and triplet states. 
The symmetrical space wave function of the singlet state may be written as 

^(1,2) = $,(1)$ B (2) + $,(2)$ B (1) + A[$,(1)$ A (2) + $ B (1)$ B (2)] (2) 

where <t> A and <t> B are the wave functions of the excitons locating at A site and B site, A is a 

parameter in the range [-1,+1]. The expectation value of energy is 

E s =< ^Ib^VZ^Ib > I < ^L^L >• The parameter X Q for the lowest energy state of the 

BE 

singlet state may be determined in a variational approach — — = . Those wave functions in 

BX 

(2) with X # X may represent excitation states. The Gaussian basis functions should be 
optimized again for the two-hole system. The antisymmetrical space wave function of the 
triplet state may be written as 

V T AB (1,2) = <D A (1)<D B (2) - <D A (2)$ B (1) (3) 

with similar definitions as in (2). The expectation value of energy is 
E T =< v I / ] iJ l//(l,2)l v I / ] iJ > / < v i / ] ij l v I , ] iJ >. The Gaussian basis set may be optimized again for 

the triplet state. The lowest expectation value of energy for both the singlet and the triplet 
states depend on the distance R AB = R A -R B as well as the dielectric constant e . The R AB 
dependence of E s and E T is shown in Figure (1) for a given value of the dielectric constant 
e = 1 . Notice that both the singlet and the triplet states are stable and the triplet state is the 
ground state for R AB > 3.5 (au.) . However, our calculation also shows that in the limit £ — > °° 
the singlet state is the ground state. The reason why the ground state depends heavily on the 
dielectric constant is that there is a strong competition between the exchange energy favoring 
the triplet state and kinetic energy favoring the singlet state. The energy gap A=E T -E S as a 
function of £ is calculated with R AB = 10 (au.) as shown in Figure (2). The gap A evolves 
on a continuum from negative value to positive value when the dielectric constant increases 
from one to infinity. For a set of parameters £ = 2.5, m* = m e , and R AB =10 (au ), the 
magnitude of gap is obtained about A = 50 mev the same magnitude of gap for high-Tc 
superconductors. Moreover, the gap is disappearing for the distance R AB > 18 (au.) as 
shown in Figure (1). The diameter of two excitons with R AB = 18 (au.) is about 
36(au.)~18A. This length scale is also the same as the coherent length in high-Tc 

superconductors. Therefore, the simple exciton model is able to give the right magnitude of 
energy gap and the coherent length, assuming that the superconducting gap be the gap 
between the triplet state and the singlet state. 

The many body wave functions for a system of N holes described by Hamiltonian H(N,M) 
with N=M>2 may be expanded in terms of the one and two-hole wave functions. For 
instance, let's consider a simple system of four holes. Four holes may form two two-exciton 
pairs, in the singlet or the triplet state. The number of configurations is 3 for two-exciton 
pairing in the singlet state, 3 in the triplet state and 6 in a mix of the singlet state and the 
triplet state. If the 3-fold degeneracy of triplet pairing is considered, there are 48 possible 
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quantum states, 48=3*1+3*9+6*3. The Hamiltonian H(4,4) may be solved in the space 
spanned by the 48 basic functions. The eigen states and eigen values may be obtained from a 
48x48 matrix of the Hamiltonian. As we have shown before, the energy level of single 
exciton split into two energy levels corresponding to the triplet and singlet state with a finite 
gap when two excitons get closer. For a system of many excitons, the two energy levels 
expands further into two energy bands due to the overlapping of wave functions. If the gap is 
larger than the band width, the ground state may be the lower band with two excitons 
pairing in the singlet state. For the simple system of 4-excitons, the wave function of the 
ground state may be expressed as a combination of the three configurations with each pair in 
singlet state. The above picture may be generalized to an N-hole system of H(N,N). The wave 
function of each configuration may take the form 

i^ {W) >=n'i* s ">n'i<^'>> w 

with \3> s u >=l v i / / '(l,2)>l^ s > and \<t> T rr >=l v I / / r r (1,2)>I^ T >. The * in the product [}' 
indicates no overlapping among the sites {//,/' J'}. The number of possible configurations is 
proportional to AM . The eigen states may be expressed generally as the combination of all the 

configurations, 1^) = V oc [u j.j.^^j vr) ) . Again for the large gap between the triplet state 

[uj'j'} 

and the singlet states, the lowest eigen states may be written as the combination of those 
configurations with those two-exciton pairs in singlet states = ^a {u] \ x ¥ {u] ) . 

[V] 

Since only those two-exciton pairs at neighbor sites are important for the ground state, the 

1*1^) may be reduced to l*P G ) = ^(x {< u>}^{<u>}^ where < IJ > indicates the neighbor sites. 

{</■/>) 

Furthermore, the coefficients cc l<u>] for those configurations {<//>} may be assumed to be 

the same, because those configurations have almost the same energy or are highly 
"degenerate". The singlet state of a two-exciton pair may be written of form 

\<t> s u > = [Cj + t c^ +c^c^ +/l(c I + t c I + L +cJ t c^)]IO), where c[ c and c^are the single-exciton 

creation operators. The ground state V¥ c ) is then expanded 

l^ G >= £ \\\<$> s 1J > = A + lB + l 2 C+-, (5) 

{</y>l {<w>) 

with A = y (ct„ c + , „ ■■■ct„ •••ct „ )I0> summing over configurations without any double 

occupation, B with only one double occupation and one vacancy, and C with only two 
occupations and two vacancies and so on. We may compare this ground state with a Bose 
condensation of mobile bosons. We use the representation given by Anderson[l]. A mobile 
boson may be written as 

b + = Y, a ( k Kr c -a = -7= L L c+ jt c+ j+ T i |0 > = -7= L c h c n |0 > ' < 6 ) 

k VM T =(„„) j V M <u> 

where nn indicates nearest neighbors. We can make a Bose condensation of N/2 bosons by 
forming 
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L T[c J ; r c + Jl t)) = A + B + C+-. (7) 



A/ 



A</2 



N/2 



<IJ> 



{<IJ>} {<IJ>} 



which is equivalent to the ground state l x P G (/l = l)>. Obviously, A = l is related to a 
vanishing Coulomb interaction in Equation (1), i.e. £ — > °° The ground state of Hamiltonian 
(1) without Coulomb interaction is a superconducting state! Usually the parameter A is 
smaller than 1 for finite Coulomb interaction. We introduce a 77 projection operation 



where n, is the occupation number operator. Acting the projection operator (8) on the state 
of Bose condensation (7) leads to a new state: P^ BS = A + (\-r])B + (\-r]) 2 C+---. The 
ground state is thus equivalent to the 77 projection of the Bose condensation, 



with 77 = 1 - A . In the strong correlation limit we have A = or 77 = 1 . The ground state 
becomes the RVB insulator defined by Anderson[l] . Moreover, *P BS is a projection of the 
BCS function *P BS = P N/2 x i i Bcs where P N/2 is a projection operator onto N/2 pairs. The 
ground state is thus an appropriate projection of the BCS function *¥ G = P ]V/2 v i / fiCS . The 77 
projection operator looks like a kind of "mobility" operator that confines the mobility of a 
Bose condensation state *¥ BS . It has to emphasize that the ground state we assume is stable 
only if the gap between the triplet and singlet states is significant. In addition, the number of 
holes on CuO may be less than the number of acceptors, N < M when the concentration of 
oxygen changes in cuprates. In this case even the RVB state might become superconducting 
[1]. The parameter A is becoming negative when two excitons are too close, corresponding to 
77 > 1 . The projection of Bose condensation becomes 

P^\m BS ) = A+\l-r)\e in B+\l-4 2 e i2n B+---. It equals to add a phase factor exp(n^) to the 
states with the n-double occupations. Further one may introduce a non-integral phase in the 
projection operator when the density of holes is very low and phase fluctuations are 
important [10]. Because the high-Tc cuprates are in metallic phase in the high doping region, 
the projected state (9) for 77 >1 might describe a metal phase. The wave function l*P G ) 
might represent either a insulator or a superconductor or a metal. It has to emphasize that 
the deduction of the equations from (5) to (9) does not depend on the dimensionality of space. 
The wave functions y ¥ c = P^ B s mav De considered as the basis functions for most of 
superconductivity. The BCS function is indeed an expansion in terms of those these basis 
functions with A = 1 and variable number of pairs. In the case of A ^ A , the Equation (9) 
represents the low-lying excitation states that are obviously still the projection of Bose 
condensation. Therefore the superconductivity is stable at finite temperature. In addition, 
due to that continuous spectrum, the specific heat should be a power function of 
temperature, rather than the exponential function for metal superconductors[ll] at low 
temperature. 

The phase diagram of H(N,N) depends on those parameters: temperature T, dielectric 
constant £ and average distance between holes R av . The last one relates to the density of 




(8) 



BS 



(9) 
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holes in terms of n = (R Q I R av ) 2 where R Q is the lattice constant for the square lattice of the 
CuO plane. In cuprates, the dielectric constant should also depend on the density of holes. 
Therefore the phase diagram depends mainly on hole density as well as temperature. In the 
low doping limit, the gap between two spin states is vanishingly small. The Hamiltonian 
H(N,N) may be mapped to a large-U Hubbard model with U »t at half-filling, 
H v = ^ J t IJ c + Ia c Ja +uy\n n n ia , where c + Ia (c Ja ) is the creation (annihilation) operator for 

Ua 1 

the bound states of excitons. The ground state of the model is an antiferromagnetic insulator 
or a RVB insulator. In intermediate doping region, the gap between the singlet and triplet 
states is significantly large. The doping dependence of the energy gap is obtained by 
A(n) = A(R) with R = R Q (l/n) 112 , R Q = 3.5 au. and e = 2.5 , as shown in Figure (3). The gap 
is obviously nonzero and positive in the region 0.05 < n < 0.45 where the ground state should 
be a superconducting state. Above Tc, the two-exciton pairs in triplet state are dominating, 
which may be responsible for the unusual normal phase properties. In even higher doping 
region, because the potential wells overlap with each other strongly, the second term of the 
Hamiltonian (1) is nearly a constant. The exciton description of holes is no longer valid. The 
Hamiltonian then reduces to a two-dimensional interacting Fermi gas that we believe should 
be metallic. In conclusion, the Hamiltonian H(N,N) has three phases: the AF or RVB 
insulator phase in low doping limit, superconducting phase in intermediate doping region; 
and metal phase in the higher doping region as shown in Figure (3). In the low doping limit 
of our phase diagram, we do not find the so-called pseudogap [12] which has a magnitude of 
1000KT. The pseudogap might be contributed by the d-holes because it has the same 
magnitude as exchange energy J. The fact that both the short-range antiferromagnetic 
order and the superconductivity coexist suggests that the d-holes might be responsible for 
the short-range antiferromagnetic order and the doping holes for the superconductivity. 

The triplet excitation states may contribute to many unusual properties of high-Tc 
superconductors. First, the probability of a singlet state being excited is proportional to 
3exp(-A / kT) instead of exp(-A / kT) as in conventional superconductors. This may explain 

the large ratio ^ r observed in cuprate superconductors. In addition, the extra number of 

states should be detected in experiments when a two-exciton pair is excited from a singlet 
state into a triplet state. Second, the triplet states may contribute to an enhanced 
susceptibility of the normal phase above Tc [13]. Furthermore, the transport properties in 
normal phase may be different from the Landau Fermi liquid, since the triplet states are 
bosons. It would be interesting to study the thermodynamics and magnetic properties of such 
a system at finite temperature. We leave them in future work. 

In conclusion, the system of doping holes in cuprate superconductors is described by a two- 
dimensional exciton model that includes the long-range Coulomb correlation. Two holes may 
form a pair in singlet state and in triplet state, the former is responsible for the 
superconducting state and the latter is responsible for the unusual normal state properties. 
The ground state is found to be the projection of a Bose condensation that may be either a 
RVB insulator, or a superconductor or a metal. The 77 projection operator introduced is 
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actually a generalization of the Gutzwiller projection operator. It is the two-dimensionality 
not the Coulomb correlation that is the driving force of superconductivity. The inclusion of 
long-range Coulomb correlation leads to a phase diagram that shows the insulator- 
superconductor transition as well as superconductor-metal phase transition. Our exciton 
model is very different from the exciton mechanism [14] that the exchange of exciton between 
two electrons is assumed to responsible for a net attractive force in analogy to the BCS 
theory. 
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Figure Captions 

Figure 1. Energy expectation value E s and E T of two-exciton pairing in singlet and triplet 
states versus the distance of the two excitons R AB . Dielectric constant is set as e = 1 . One 
energy unit equals to 27.2 ev. 

Figure 2. Energy gap A = (E T - E s ) between the triplet and the singlet states versus 
dielectric constants . The distance between two excitons is set as R AB = 10 au.. 

Figure 3. Phase diagram of the Hamiltonian (1) and the doping dependence of gap for a 
given dielectric constant e = 2.5 and lattice constant R Q = 3.5 au. . 
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Table 1. The results of numerical calculation for a single exciton and a Hydrogen atom in the 
variational approach. 



Wave function c, (unrenormalized), optimized Gaussians a, , expectation value of energy E 
and mean radius < r > (in atomic unit) for exciton and Hydrogen atom are listed. 





a, 


E(Hartrees) 


c t 


a, 


E(Hartrees) 


(exciton) 


(exciton) 


1 <r> (au.) 


(Hydrogen) 


(Hydrogen) 


1 <r> (au.) 


c l = 1.0 


a! = 0.0198 


-0.159 


cj = 1.0 


a t =0.159 


- 0.497 


c 2 = 1.04 


a 2 = 0.0582 


/ 7.228 


c 2 = 1.843 


a 2 = 0.775 


/l 


c 3 = 0.0056 


a } = 3.3281 




c 3 = 1.168 


a 3 = 7.54 
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Figure 2 




Figure 3 



Phase Diagram 
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Density of holes n=(R o /R) , R o =3.5au., t =2.5 
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